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Transport coefficients of the second order hydrodynamics 
and the applicabihty of hydrodynamic model 

Shin MuROYA 

Dept. of CM. Matsumoto Univ. Matsumoto 390-1295, Japan 

Based on the Nakajima-Zubarev type nonequilibrium density operator, we derive micro- 
scopic formulae of the transport coefficients in the second order hydrodynamics. 

§1. Introduction 



Hydrodynamical model is one of the widely applied phenomenological models 
for the relativistic heavy ion collisions. According to the recent detailed analyses, 
g3 \ weak but non-vanishing viscosity is significant for the qualitative discussion of V2 

Q ■ at RHICP Because Navier-Stokes equation is a parabolic type equation, naive rel- 

ativistic extension is not consistent with causality and numerical solution becomes 
unstableP The relativistic causal hydrodynamics had been introduced phenomeno- 
logically by Israel and Stewart more than thirty years ago.'^'SI' In these decades, 
many proposal appear for the appropriate equations of the relativistic causal hydro- 
dynamicsP®''^ Most of the works seem to belong semi-classical phenomenological 
CJ ' approach based on Boltzmann equation. In this paper we discuss a relativistic causal 

2 ■ hydrodynamics as the second order hydrodynamics based on a nonequilibrium den- 

sity operator method-'^' 

Hydrodynamical model is a phenomenological model of the macroscopic point 
<^ . of view. By coarse graining, system is constructed as a patchwork of corresponding 

Q^ \ statistical systems in equilibrium. The system at the space-time point is micro- 

■ scopically large but macroscopically one point and moving with four velocity . 

Thermodynamical quantities and coefficients of hydrodynamical equation should be 



the local quantities as the function of position through the local thermodynamical 
parameters such as temperature T[x) and chemical potential In principle, we 

can calculate them by using statistical mechanics. The aim of this paper is to de- 
scribe all coefficients in the second order hydrodynamics as expectation values of the 
local equilibrium state. 



■ §2. Notations 

Hydrodynamical equation is composed of energy-momentum conservation low, 
d^Tpfj = 0, and charge conservation low, d^Jp = 0,. Flow of the fluid is described by 
a four velocity C/'* which is a normalized time-like vector, U^Up = 1, and by using 

space- like projection operator is defined as, A^^'^ = g^" — U^V^ . 

We adopt Landau-Lifshitz frame (E- frame); is eigenvector of the energy 
momentum tensor,!^ T^^^Uy. = ell^. Pressure P is defined by the trace of spatial 
part of the energy-momentum tensor, P = ^Z\^^T^'^. Stress-shear tensor is given 
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as the spatial traceless part, 

Charge density is given by n = J^Ufj, as usual, however, is proportional to 
the energy flow on the E-frame, charge current contains P'' which is perpendicular 
to 

T stands for time on the comoving flame of the fluid element, r = x^U^, and D = 
U^d^ stands for the derivative by r, respectively. 

§3. Nonequilibrium density operator 

Let us focus our discussion to the only small deviation from thermal equilib- 
rium state. Our aim is to derive equations which enable us to discuss second order 
deviation from the thermal equilibrium state. 

Nonequilibrium density operator is given as,'Sl'''2} 

p = Q~^exp 



dx^ {/3(a;, t)U''{x, t) ^(a;, t) -^{x, t)fi{x, t)J^{x, t)} 
lim^y'd^aj'e-^^''^^-^')'' ^f^^,{x,' t')^'^'{P{x' ,1')^ {x\t')) 

C/^(x - x')^ > 0. (3-1) 



C 

+ Ux',t')d"'{fi{x\t')p{x',t'))^ 



The first term of the right hand side corresponds to the local equilibrium density 
operator, pi = Qj'^ex.p{J dxf^{l3{x, t)U'^{x, t)T^y{x, t) — I3{x,t)fi{x, t)J^{x, t)}), with 
Ql being normalization, dx^ is a space like hypersurface which locally corresponds 
to a comoving flame of the fluid element. Though out this paper, the limit, ( — > 0+, 
should be taken at the final stage of the calculation. The second term stands for the 
effect of thermodynamical force. Expectation value of any operator 0{x) is given by 
(d(x)) = tr{p6(x)}. 

Following the standard procedure of linear response theory, we expand p around 
local equilibrium operator, pi, up to the linear term of thermodynamical force. Ex- 
pectation value, (0(x)), is given as a local equilibrium expectation of canonical 
commutation with the force term, 



= {d{x))i + ^lirn^y fiVe-^^''(^-^')'' (d{x),fpAx')d"' ))) 

+ lim^ / d^a^'e-^^''^^-^')^ (d{x), Jp{x')d"' {p{x')p{x'))^ , (3-2) 



with {0,B) being canonical commutation relation 



iO,B) = {O [\xe'^^Be~^^)i - {0)i{B)i. 
Jo 
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§4. Derivative expansion with thermodynamical parameters 

The change of thermodynamical parameters are slow enough and almost con- 
stant during the microscopic relaxation length which is determined by the canon- 
ical correlation, we may adopt Taylor expansion of the thermodynamical force, 
d'P {l3{x')U" {x')) and d'^ {/3{x')fj,{x')) at around position x of the operator 0{x). 
On the local comoving flame of the fluid element at x, we can replace the local 
equilibrium distribution with the corresponding distribution in equilibrium. 

Then we can obtain expectation values of T^'^ and in the power series of 
gradient of thermodynamical parameters, 

{f^"'{x)) = {f^''{x))i + j d^x' f (it'e-^(*-*')(f^^(x), {pU") 

J J — oo 

+ I d^x' f dt'e-^(*-*')(f^^(x), Jp(x'))a''(/3/u) 

J J —OO 

+ [ d'x' f dt'e-^(*-*')(f^^(x), [xx' - xx)fp,{x'))d^dP {f3Un 

J J —-CO 

+ / d'x' f dt'e-^(*-*')(f ^^(x), [xx' - xx)Jp{x'))d^dP (/3/i) , (4-1) 

J J — oo 

= {>{x))i + / d^x' f (it'e-'^(*-*')(>(x), Jp{x'))dP iPn) 

J J — OO 

d^x' f dt'e~^^'"''\j^{x),fp„{x'))dP (/?[/") 

J —oo 

d^x' r di'e-'^(*-*')(>(x), {xx' - xx)Jp{x'))d^dP (/3^) 

J — oo 

+ / d'x' f dt'e-^'^'-''\jP{x), {xx' - xx)fp,{x'))d^dP {pun .(4-2) 



-oo 



The first terms of the right hand side of (|4-1|) and (|4-2p are the zeroth order terms 
which correspond to the perfect fluid part, {TP'^{x))i = (e -|- P)UPU'^ — PgP'^ , and 
{Jp{x))i = nJP. The second terms of the right hand side provide us Navier-Stokes 
equation. Well known Kubo-formula for viscosity and heat conductivity are obtained 
as the coefficients of the first order gradient of thermodynamical parameters .l^'l^ 

In both equations, the third terms vanish by virtue of Curie's theorem, because 
of isotropic property of the local equilibrium state. The forth and the fifth terms 
stand for the second order deviations. For example, the forth term in right hand 
side of (|4-2p can be rewritten as, 

d^x' t (it'e-^(*-*') {xx' - xx)Jp{x')) d^dP {pfi) 

-oo 

d^x' I dt'e-^(*-*')(>(x), {t' - t)Jp{x'))DdP {Pfi) 
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+ j (fx' f dt'e-^^*-*'\j^{x),A''^{xx' -xx)Jp{x'))A„pd^dP{l3fi). (4-3) 

J J — oo 

Comparing (|4-2p and (2.38b) in Ref. 4, the first term of tlie right hand side of 
(|4-3p corresponds to the time-derivative term of heat current and the coefficient 
corresponds to K?T^n? Pi/ {e + P)^ in Israel- Stewart equation-Sl On the other hand, 
the second term in ()4-3p must vanish by virtue of Curie's theorem. 



§5. Concluding Remarks 



All coefficients in the derivative expansion in (j4-ip and ()4-2p are expressed by 
the canonical correlations of T^^ ^ and x'^ — x^. Coefficients for the second order 
terms are current-weighted relaxation time, 

/ d^x' f dt'e~^(*-*')(>(x), {t' - t)Jp{x')), 

J J — oo 

and current-weighted correlation length, 

/ d^x' f dt'e-^'^^~^'\f^'^{x),A''^{x)! -xx)Jp{x')). 

J J — oo 

These quantities can be evaluated by using hadro-molecular simulatioiP2)>[I31ior Lat- 
tice field theory. Comparing these quantities to the scale of the change of ther- 
modynamical parameters in the hydrodynamical solution, we can investigate the 
self-consistent check of the hydrodynamical picture.l^ 
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